从光闪烁频谱反演大气湍流谱:原理与数值模拟

建立了一种数值算法,从光波在湍流大气中传播的闪烁时间频谱的高频特征反演大气光学湍流谱.此算法具有两点优越性:（1）反演出的湍流谱直接表征大气折射率的起伏而非大气温度场的起伏;（2）可以获得高空间波数、大部分位于耗散区的湍流谱,而用其它方法不易获得此区域的湍流谱.对Hill湍流谱以及一种复杂的湍流谱进行了大量数值模拟来检验此种算法,结果表明:如果知道大气风速与湍流内尺度就可以相当有效地获得湍流谱.     

MACRO AND MICROSCOPIC PARAMETER MAPPING RELATIONSHIP FOR RED SANDSTONE 1)

颗粒流离散元软件PFC2D在岩土工程中的应用十分广泛,存在的主要问题是如何标定其细观参数,目前大多使用"试凑法",此方法的缺点在于工作量大、效率低等,本文提出以反向传播算法(back propagation,BP)神经网络的方式代替此方法,利用PFC2D 内置的FISH 以及Python 语言对其进行二次开发,使之自动运行并获取40组宏细观参数样本。结果表明：BP神经网络可以快速准确地建立宏细观参数映射关系,误差均保持在0.01 以内,模拟得到的应力--应变曲线与室内试验曲线高度吻合,且无需大量的数据样本便可创建网络,效率较高;另外,经验证本文选用的平直节理模型,可以有效地解决平行粘结模型UCS/TS值偏小问题,确定平直节理模型可以更好地模拟岩石。     

动态裂纹扩展中的分形效应

假设裂纹顶端沿着分形轨迹运动，建立了裂纹扩展的分形弯折（ｋｉｎｋｉｎｇ）模型来描述裂纹的动态扩展。根据这个模型，我们推导了分形裂纹扩展对劝态应力强度和裂纹速度的影响．动态应力强度因子与表观应力强度因子之比Ｋ（ｌ（ｔ），Ｖ）／Ｋ（Ｌ（ｔ），Ｏ）是表观裂纹速度Ｖ＿Ｏ，材料微结构参数（ｄ／Δａ），分维Ｄ和裂纹扩展路径的弯折角θ的函数。本文研究结果表明：在分形裂纹扩展中，表观（或量测）的裂纹速度Ｖ＿Ｏ很难接近Ｒａｙｌｅｉｇｈ波速Ｃ＿ｒ．动态断裂实验中Ｖ＿Ｏ明显低于Ｃ＿ｒ的原因可能是分形裂纹扩展效应所致。材料的微结构，裂纹扩展路径的分维和弯折角均很强地影响动态应力强度因子和裂纹扩展速度。     

一种岩石损伤本构模型在地下强爆炸中的应用

为了描述地下爆炸波的传播,建立了一种考虑塑性硬化、剪胀和损伤软化效应的岩石本构模型,并应用于地下强爆炸自由场的数值计算,获得的速度和位移波形以及峰值速度衰减曲线等与国外地下试验测量数据和同类计算相比十分接近,从而验证了本构模型的有效性。     

Application of scaled boundary finite element method in static and dynamic fracture problems

The scaled boundary finite element method (SBFEM) is a recently developed numerical method combining advantages of both finite element methods (FEM) and boundary element methods (BEM) and with its own special features as well. One of the most prominent advantages is its capability of calculating stress intensity factors (SIFs) directly from the stress solutions whose singularities at crack tips are analytically represented. This advantage is taken in this study to model static and dynamic fracture problems. For static problems, a remeshing algorithm as simple as used in the BEM is developed while retaining the generality and flexibility of the FEM. Fully-automatic modelling of the mixed-mode crack propagation is then realised by combining the remeshing algorithm with a propagation criterion. For dynamic fracture problems, a newly developed series-increasing solution to the SBFEM governing equations in the frequency domain is applied to calculate dynamic SIFs. Three plane problems are modelled. The numerical results show that the SBFEM can accurately predict static and dynamic SIFs, cracking paths and load-displacement curves, using only a fraction of degrees of freedom generally needed by the traditional finite element methods.The project supported by the National Natural Science Foundation of China (50579081) and the Australian Research Council (DP0452681)The English text was polished by Keren Wang.     

阶梯压电层合梁的波动动力学特性

采用行波理论系统地研究了压电阶梯梁的自由振动分析以及强迫响应的分析方法. 基于分布 参数理论研究了压电阶梯梁的波传播特性，忽略柔性梁横向剪切和转动惯量的影响，给出了 梁的轴向和横向的简谐波解. 将压电阶梯梁离散化为单元，考虑压电片的刚度和质量的影响， 建立了节点散射模型. 应用位移连续和力平衡条件，推导了节点的波反射和波传递矩阵，在 此基础上，引入波循环矩阵的概念，给出波循环矩阵、波传递系数矩阵的确定方法. 应用波 循环矩阵可以有效地计算结构的固有频率. 另外，应用波传递系数研究了压电陶瓷作动器位 置对其驱动能力的影响. 得出两个主要结论：1)作动器靠近悬臂梁固定端将有较强的驱动 能力，悬臂梁边界反射行波产生弯曲消失波有利于增大压电波的模态传递系数；2)模态传 递系数与固有频率的灵敏度密切相关，波传递系数越大, 对应该处固有频率变化灵敏度越大. 另外，数值算例表明了行波方法比有限元方法具有更高的计算精度.     

Theoretical and experimental analysis of longitudinal wave propagation in cylindrical viscoelastic rods

Wave propagation in viscoelastic rods is encountered in many applications including studies of impact and fracture under high strain rates and characterization of the dynamic behavior of viscoelastic materials. For viscoelastic materials, both material and geometric dispersion are possible when the diameter of the rod is of the same order as the wavelength. In this work, we simplify the Pochhammer frequency equation for low and intermediate loss viscoelastic materials and formulate corrections for geometric dispersion for both the phase velocity and attenuation. The formulation is then experimentally verified with measurements of the phase velocity and attenuation in commercial polymethylmethacrylate rods that are 12 and in diameter. Without correcting for geometric dispersion, the usable frequency range for determining the phase velocity and attenuation for the rod is about , and about for the rod. Using the correction procedure developed here, it was possible to accurately determine the phase velocity and attenuation up to frequencies exceeding for the rod and for the rod. These corrections are applicable to many polymers and other viscoelastic materials. From thereon, the viscoelastic properties of the material can be determined over a wide range of frequencies.     

增韧环氧树脂的动态裂纹扩展研究

本文主要进行了环氧及增韧环氧树脂的断裂韧性及裂纹快速扩展的试验研究。试验过程中采用了ＧＬＣ－１型高速裂纹扩展测试仪来测试裂纹的扩展速度,得到在裂纹扩展过程中裂纹扩展速度曲线。本文结合不同的计算公式及有限元分析方法,讨论了各个确定断裂韧性公式的准确程度,发现传统的静态断裂韧性的分析方法所得到的结果偏大,有一定的危险性,建议使用试验与数值计算相结合的方法;同时还发现增韧不仅可以提高材料的静态和动态断裂性能,而且在裂纹扩展过程中可以起到减缓裂纹扩展的作用     

Seismic gyrotropy as a result of rock-microstructure dissymmetry

Seismic investigations show that rocks exhibit not only anisotropic properties but also gyrotropic ones. Recently, the phenomenological theory has been developed to describe the phenomenon ‘seismic gyrotropy'. The paper gives a possible explanation of the gyrotropic properties of rocks on a microlevel. A gyrotropic model of grainy rock with dissymetry in its microstructure is presented. The model is built on a principle ‘azimuthal turn plus translation'. The mathematical and physical formulations of the problem are discussed. An algorithm is constructed to find the displacements, their first and second derivatives with respect to spatial coordinates as well as stresses inside a grain, an element of the model. Mathematical modelling of elastic-wave propagation in the model considered enabled to demonstrate that the model possesses gyrotropic properties, i.e., it rotates on a polarization vector. Gyration constants are determined.     

Dynamic compression-shear response of brittle materials with specimen recovery

A new configuration for compression-shear soft-recovery experiments is presented. This technique is used to investigate various failure mechanisms during dynamic multiaxial loading of an Al₂O₃/SiC nanocomposite and TiB₂. Velocity profiles of the target surface are measured with a variable sensitivity displacement interferometer, yielding normal and transverse velocity-time histories. A dynamic shear stress of approximately 280 MPa is obtained, in the Al₂O₃/SiC nanocomposite, for an imposed axial stress of about 3.45 GPa on a 540 m thick sample. This dynamic shear stress is well below the value predicted by elastic wave propagation theory. This could be the result of stress-induced damage and inelasticity in the bulk of the sample or inelasticity on the sample surface due to frictional sliding. To gain further insight into the possible failure mechanisms, an investigation of compression-shear recovery techniques, with simultaneous trapping of longitudinal and lateral release waves, is conducted.